Indications of Unconventional Superconductivity 
in Doped and Undoped Triangular Antiferromagnets 
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The possibility of superconductivity in doped and undoped triangular antiferromagnets is dis- 
cussed. Using the Bethe-Salpeter (B-S) equation, it is shown that the exchange of RPA paramagnons 
on a triangular lattice Hubbard model leads to strong pairing correlations at and near half-filling. 
The dominant states for this system correspond to d-wave singlet (even-frequency) and s-wave triplet 
(odd-frequency) pairing. Analytical techniques applied to the hole-doped t-J model yield similar 
results. A ti — t2 Hubbard model interpolating between square (ti = 0) and triangular (ti = £2) 
lattices has a tendency to only d-wave singlet pairing for t\/t2 < 0.8. Experimental consequences 
for organic compounds ac-(BEDT-TTF)2X are discussed. 



Since the discovery of high-temperature supercon- 
ductivity, two-dimensional (2D) doped antiferromagnets 
(AF) have been studied intensively. Experimental re- 
sults and a large body of theoretical work suggest that 
the pairing state is highly anisotropic and probably a 
d 3 .2„j / 2-singlet. In addition, materials with 2D spin ar- 
rangements on non-square lattices have also been syn- 
thetized. For example, experiments suggest, the realiza- 
tion of a triangular spin-i AF in NaTiO^, as well as 
in surface structured^ such as K/Si(lll):B. Delafossite 
cuprates RCu0 2 +,5, with R a rare-earth element, have 
Cu ions sitting on a triangular latticeS. Otheij-compounds 
with this geometry have also been preparedcl. Recently, 
theoretical arguments in the context of organic supercon- 
ductors suggest that k-(BEDT-TTF) 2 X, where X is an 
ion, is described by the half-filled Hubbard model on 
anisotropic triangular lattice with one electron per sit 
Metallicity is obtained by working below the critical cou- 
pling that leads to an insulating phase. 

The presence of unconventional superconductivity on 
a square-lattice AF suggests that similar unusual phe- 
nomena may occur in other 2D AF arrangements, such 
as a triangular lattice. Although most of the materials 
described above do not contain a finite concentration of 
mobile carriers, it is important to analyze on theoretical 
grounds the properties of holes in this environment and 
specially the characteristics of a possible superconduct- 
ing state. This information may induce further exper- 
imental work in non-square lattice doped spin systems. 
Precisely in this paper the issue of pairing in triangular 
and anisotropic triangular AF is addressed. 

The model to be discussed in the first part of this 
paper contains a nearest-neighbor hopping term i 2 on 
a 2D square lattice and an additional hopping term t\ 
along one diagonal of each plaquette, so that t± = 
corresponds to a square lattice (bandwidth 8i 2 ) whereas 
ti = i 2 is an isotropic triangular lattice (bandwidth 9t 2 ). 
For details see Ref.u. The interaction will be modeled by 
an on-site Hubbard repulsion Un^rii^. 

To study the effective interaction between electrons 
arising from the exchange of paramagnons the RPA ap- 



proximation is used, as done n before for the analogous 
problem on the square lattices where further worku us- 
ing more sophisticated techniques showed that indeed in- 
dications of d-wave pairing are present in this context. 
Then, the B-S approach is powerful enough to capture 
the essence of the pairing tendencies in Hubbard-like 
models. The effective pairing potentials for the singlet 
and the triplet channels are 
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with p = (p, iu) n ) and uo n being a fermionic Matsubara 
frequency. Xo(q>w) denotes the spin susceptibility: 
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where e p are the non-interacting electron energies. 
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FIG. 1. Effective singlet (V s ) and triplet (Vt) interaction 
in real-space for U/t2 = 4, ti = ti, uo = and a density of 
(n) ~ 0.9 (results at half-filling are similar). The numbers are 
normalized such that the on-site interaction has unit magni- 
tude. (The reference site is marked by a rectangle.) The ratio 
of on-site interactions is V 3 (0) /Vt(0) ~ —3. With decreasing 
band filling the range and magnitude of the interactions de- 



These interactions Fourier-transformed into real-space 
at ti — i 2 (triangular lattice) and in the static limit 
{u> = 0) are shown in Fig. |l]. The density is chosen to 
be slightly below half filling, i.e., (n) « 0.9. The singlet 
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pairing interaction is repulsive on-site, but attractive for 
pairs of electrons at nearest-neighbor sites. The triplet 
interaction Vt behaves qualitatively similar to V s except 
for the difference in sign, i.e., it is repulsive for nearest- 
neighbor pairs. Note that the triangular system does not 
obey perfect nesting at half-filling. Thus, the tendency 
to form a spin-density wave state is not as strong as on 
the square lattice Hubbard model. However, for large U 
and half-filling the model maps onto the S = | triangular 



-believed to have an ordered 120° magnetic 



AFO which is 
ground states. 

To investigate two-particle bound states it is conve- 
nient to classify them according to symmetry proper- 
ties in momentum space. The point symmetry group 
of the isotropic triangular lattice is Cq which has 6 irre- 
ducible representations (s, p x , p y , d x 2_ y z, d xy , f x 3-3 X y 2 )- 
The results show that the triangular system has possible 
instabilities in the s-, p x -, d xy - and / 2 ,3_ 3:I ,j ; 2-channels 
illustrated in Fig. || using the Brillouin zone (BZ) 
of the triangular system. There are 3 equivalent d xy - 
solutions (related by 60° rotation) , one of them continu- 
ously evolves into the well-known square-lattice d x i_ y i- 
symmetry when reducing t\ and going to the square BZ. 
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FIG. 2. Momentum dependence of the pair wavefunctions 
corresponding to s, p x , d xy , and f x 3-: ixy 2 symmetries. The 
graphs show the sign and the zeros of the pair wavefunction 
at selected points of the BZ. 

To search for possible pairing instabilities in the prob- 
lem let us here investigate the leading eigenvalues A and 
eigenfunctions <fi(p) of the Bethe-Salpeter (B-S) equation 
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Self-energy contributions are neglected, and the single- 
particle Green's function is G(p) — (iuj n — e P ) _1 . When 
the largest eigenvalue A reaches 1, a superconducting in- 
stability to a state with a pair wavefunction <f)(p) occurs. 

The B-S equation can have singlet and triplet so- 
lutions which correspond to pair wavefunctions with 
even or odd parity when transforming from (p, iu) n ) to 
(— p, —iui n ). The usual s- and d-wave singlet solutions 
are even both in frequency and momentum whereas the 
usual triplet solutions (p and /) are even in frequency 
and odd in momentum. As discussed by Berenzinskii 
and by Balatsky and AbrahamsEHl also odd-frequency 
pair wavefunctions are permitted by symmetry. Then, 
one can obtain even-momentum triplets ((/)(— p,itu n ) — 
4>(p, iuJn) — — 0(p, — iijJn)) and odd-momentum singlets 
(4>(—p,icj n ) = -<p(p,iu) n ) = 4>(p,-iw n )). In a model 



with strong on-site repulsion there are two ways for the 
pairs to "avoid" the repulsion. One is by construct- 
ing a pair wavefunction with zero on-site amplitude, i.e. 
with non-zero angular momentum. The second way is by 
having an odd-uj dependence which implies a vanishing 
equal-time pair amplitude (and spontaneous breaking of 
time reversal symmetry). 

The present investigation shows that the dominant 
pairing channels in the triangular lattice model treated 
with the B-S equation are d-wave singlet (even in fre- 
quency) at intermediate U/t2 switching to s-wave triplet 
(odd in frequency) at large U/t2, although a small de- 
viation from this geometry using the t\ — £2 model es- 
tablishes a d-wave state as the only dominant one. (For 
small XJ jt% < 2.5 and t\ = £2 the transition tempera- 
ture is smaller than 10~ 3 £2-) This curious competition 
for the triangular lattice can be understood intuitively 
as follows, starting with the even-w case first: The sin- 
glet potential favors pairing with zero on-site amplitude 
due to the Hubbard repulsion. This effect suppresses s- 
wave singlets, as is the case in the square lattice near 
half-fillingtx In the family of even-w solutions the next 
possibility is a d-wave singlet. However, in contrast to 
the square lattice where d^^^-pairing is favoredEl, the 
tendency towards d-wave solutions is not as strong in 
the triangular system. This can be understood from the 
B-S equation: The susceptibility xo an d the effective po- 
tential are strongly peaked at the ordering wavevector 
Q = (~tt, 0) for the triangular system. The kernel of 
the B-S equation for T — > basically couples values of 
4>(p) at Fermi surface momenta which differ by Q. The 
particular value of Q for the triangular system weakens 
d x 2_ y 2- and d xy -wa,ve singlet solutions. How about p and 
/ symmetries in the even-w sector? Here the effective 
triplet potential is repulsive for nearest neighbors, and it 
tends to suppress all triplet pairing with zero amplitude 
on-site. This reasoning lead us to believe that odd-w 
pairing may be relevant for the triangular lattice. In the 
triplet channel, the s-wave solution can take advantage of 
the on-site interaction which remains attractive at non- 
equal times, thus avoiding the suppression by the Pauli 
principle. Note that in previous B-S studies on the square 
lattice it was observed that odd-w solutions actuallypGpm- 
pete in strength with the even-w d-wave possibility!!! 

To verify this intuitive picture, the B-S equation was 
solved on a finite lattice (up to 32 2 sites) with a cutoff on 
the Matsubara frequencies of 8 £2- The eigenvalues reach 
unity only at and close to half-filling. The transition tem- 
perature is determined by the temperature where this 
largest eigenvalue reaches unity. In Fig. || the leading 
eigenvalues are shown for the triangular lattice (£1 = £2) 
at half-filling, U/t = 3.5, vs temperature. In this regime 
the dominant eigenvalue corresponds to d-wave singlet 
pairing. Small deviations from the perfect triangular lat- 
tice using the £1 — £2 model also give a stable d-wave 
singlet solution at a temperature T = 0.08 £2 (Fig. pi). 
This regime is likely of relevance for organic supercon- 
ductors^. However, it is interesting to observe that for 
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U/t > 4 and half-filling the d-wave even-w and the three 
odd-w solutions strongly compete in the considered tem- 
perature range O.Of < T/t 2 < O.f (Fig. |b). Actually, 
the eigenvalue corresponding to the s-wave pairing insta- 
bility is the largest in this regime, and thus odd-uj s-wave 
triplets are favored, as expected from the previous dis- 
cussion. Since self-energy contributions are neglected the 
critical temperature for the pairing instability, which is 
around 0.02 1 2 , should only be considered as a rough es- 
timate—Compared with the results obtained on a square 
latticeO (~ 0.1^), the triangular lattice seems to have 
a critical temperature substantially smaller. 
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FIG. 3. Bethe-Salpeter eigenvalues vs temperature for 
the isotropic triangular lattice, U = 3.5 t2, and half-filling 
(d-wave singlets solid, s-wave triplets dash-dot, p-wave sin- 
glets dashed, /-wave singlets dotted). 
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FIG. 4. Bethe-Salpeter eigenvalues vs. band filling for 
U = 4t 2 . a) txjti = 0.8 and T = 0.08i 2 - b) ti/h = 1.0 
(triangular lattice) and T — 0.02 t2- The inset in a) shows 
the symmetry of the d-wave solution in the triangular lattice 
BZ which becomes d x 2_ y 2 for the square lattice case. 

Thus far, the calculations have addressed the weak 
coupling limit. In the second part of this paper an 
isotropic triangular AF doped with holes described by 
the t-J model will be discussed. At half-hlling the t-J 
model reduces to the triangular Heisenberg antiferromag- 
nct (THAF)lj with spins i. It is by now widely be- 
lieved,that its ground state possesses magnetic long-range 
orderO which can be described as 120°-order, namely the 
ground state of the classical spin model with the addition 
of quantum fluctuations. 

The one- hole problem in the triangular t-J model was 
studied analytically befoffi^using the self-consistent Born 
approximation (SCBAjtlia Since the SCBA neglects 



closed loops in the hole motion (which are more impor- 
tant here than for the square lattice since only three hop- 
ping steps are necessary for one loop) in this paper an 
alternative analytical approach based on the picture of 
magnetic polarons will be used. The spin excitations are 
described by a set of string operators, locally distorting 
the background spin conhguratiorJi-ali3 which is assumed 
to be long-range ordered. The one-hole spectral func- 
tion can be evaluated using the MarhZwanzig projection 
technique. For details see Refs.E£EZI. Fig. |5|a contains 
the calculated one-hole dispersion £p 0l ° for different cou- 
plings J/t. In Fig. ||b the one-hole spectral function at 
J/t = 0.4 and the momentum corresponding to the band 
minimum p = (|tt, 0) is shown. A sharp peak at the 
bottom of the spectrum suggests that a description in 
terms of quasiparticles (QP) is appropriate. Note that 
the one-hole dispersion shown in Fig. ||a is in qualita- 
tive agreement |sdth exact diagonalization (ED) data on 
a 12-site clusteniias well as with SCBA and 21-site ED 
results from Ref£3. (A detailed discussion will be given 
elsewhercO.) 
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FIG. 5. a) One-hole dispersion calculated from the po- 
laron picture at J/t = 0.1 (dashed) and 0.4 (solid). The 
energy zero level has been set at the center of mass of the 
band. The inset shows the path in momentum space corre- 
sponding to the horizontal axis, b) One-hole spectral function 
at J/t — 0.4 and momentum (|7r, 0) which is the minimum of 
the quasiparticle (QP) band. A sharp QP peak is observed. 

Now let us consider a fermionic model constructed 
from the dressed one-hole properties and including an 
interaction between the QP. The use of a dressed QP 
dispersion implicitly assumes that in the normal state 
short-range antiferromagnetism exists, namely the AF 
correlation length is not negligible. This is a standard 
assumption in AF scenarios for the cuprates and can be 
checked experimentally using, e.g., NMR techniques. For 
simplicity, the QP interaction is taken from the RPA cal- 
culation described in the first part of the paper. For sin- 
glet pairing it basically consists of an on-site repulsion, a 
nearest-neighbor attraction, and a NNN repulsion. Note 
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here that the w-dependence of the interaction is impor- 
tant for possible odd-w pairing 1 J'revious literature in the 
context of the square lattice LJfiJ has shown that in both 
weak and strong coupling the most important part of the 
interaction for the study of pairing is the NN attraction 
(singlet channel). Thus, for the qualitative purposes of 
the present strong coupling study it is reasonable to as- 
sume that a similar result holds in triangular lattices. 

To investigate the existence of pairs within this model 
once again the B-S equation Eq.(|) will be solved. Now 
the single-particle Green's function is determined by the 
one-hole QP dispersion from the strong-coupling calcula- 
tion: G(p) = (iUn-e^)- 1 . For the calculation the RPA 
interaction V(p,p') corresponding to an electronic den- 
sity (n) = 0.9 will be used. In Fig. || the leading eigen- 
values are plotted as a function of temperature for the 
different pairing symmetries, working with a low-band 
hole filling (10%). The results qualitatively agree with 
those obtained in the Hubbard model. 



the isotropic ones at intermediate U/t. However, at large 
U jt pairing in the s-wavc triplet odd-frequency channel is 
also possible. Hopefully the present calculations will mo- 
tivate experimentalists to study in more detail antiferro^. 
magnetic triangular lattice compounds. Recent resultsEl 
suggest that the physics described here may already be 
realized in organic superconductors_If ARPES experi- 
ments were possible in this contextEll, nodes should be 
observed for a d-wave singlet state.. Specific heat and 
penetration depth measurementsoEj already suggest the 
existence of gapless excitations in these materials. 
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FIG. 6. Temperature dependence of the Bethe-Salpeter 
eigenvalues for s-wave triplets, and p-wave, d-wave singlets, 
and /-wave singlets, now calculated from the strong-coupling 
dispersion and the RPA interaction. 

Odd-frequency pairing such as the one described here 
has experimental interesting consequences, as discussed 
in Ref.Ej. First, there is no gap in the excitation spec- 
trum at any k. Second, the magnitude of the Hebel- 
Slichter peak in the NMR relaxation rate is strongly re- 
duced compared to the BCS case. Other nontrivial con- 
sequences of odd-frequency pairing include the tunneling 
between an odd-frequency and a BCS superconductor. 

However, regarding the organic superconductors k- 
(BEDT-TTF) 2 X, our results suggest that their low- 
temperature properties will likely correspond to a d-wave 
singlet state being the analog of the widely discussed 
d x 2_ y 2 -state of square lattices which arises from antifcr- 
romagnetic fluctuations. 

Summarizing, in this paper the possibility of su- 
perconductivity in doped and undoped triangular and 
anisotropic antiferromagnets has been investigated. The 
pairing interaction is assumed to arise from the ex- 
change of paramagnons. Both weak-coupling and strong- 
coupling approaches suggest that the system presents a 
pairing instability in an unconventional channel. A d- 
wave singlet superconducting state similar to the square 
lattice dominates for anisotropic lattices, and even for 
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